Enhanced Transmission of Light and Matter through Nanoapertures without 

Assistance of Surface Waves 
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Subwavelength aperture arrays in thin metal films enable enhanced transmission of light and 
matter waves [for example, see T.W. Ebbesen et al., Nature (London) 391, 667 (1998) and E. Moreno 
et al., Phys. Rev. Lett. 95, 170406 (2005)]. The phenomenon relies on resonant excitation of the 
surface electron or matter waves. We show another mechanism that provides a great transmission 
enhancement not by coupling to the surface waves but by the interference of diffracted evanescent 
waves in the far-field zone. Verification of the mechanism is presented by comparison with recently 
published data. 

PACS numbers: 42.25.Bs, 42.25.Fx, 42.79.Ag, 42.79.Dj 
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The scattering of waves by apertures is one of the ba- 
sic phenomena in the wave physics. The most remarkable 
feature of the light scattering by subwavelength apertures 
in a metal screen is enhancement of the light by excitation 
of electron waves in the metal. Since the observation of 
enhanced transmission of light through a 2D array of sub- 
wavelength metal nanoholes , the phenomenon attracts 
increasing interest of researchers because of its poten- 
tial for applications in nanooptics and nanophotonics . 
Recently, the enhanced transmission through subwave- 
length apertures was predicted also for matter waves [sj . 
The enhancement of light is a process that can include 
resonant excitation and interference of surface plasmons 
H H H H , Fabry-Perot-like intraslit modes 0, H H E3 > 
and evanescent electromagnetic waves at the metal sur- 
face ■ In the case of thin screens whose thickness are 
too small to support the intraslit resonance, the extraor- 
dinary transmission is caused by the excitation of surface 
plasmons or their matter-wave analog, surface matter 
waves I^^IMI^' In this Letter, we show another mech- 
anism that provides a great transmission enhancement 
not by coupling to the surface electron or matter waves 
but by the interference of diffracted evanescent waves in 
the far-field zone. 

The transmission enhancement without assistant of 
surface waves can be explained in terms of the follow- 
ing theoretical formulation. We first consider the trans- 
mission of light through a structure that is similar but 
simpler than an array of holes, namely an array of par- 
allel subwavelength slits. The structure consists of a 
thin metal screen with the slits separated by many wave- 
length. To exclude the plasmons from our model, the 
metal is considered to be a perfect conductor. Such a 
metal is described by the classic Drude model for which 
the plasmon frequency tends towards infinity. Owing to 
the great distance between the slits, the electromagnetic 
field at one slit is assumed to be independent from other 
slits. The transmission of the slit array is determined 
by calculating the light power in the far-field diffraction 
zone. The waves diffracted by each of the independent 



slits are found by using the Neerhoff and Mur approach, 
which uses a Green's function formalism for a rigorous 
numerical solution of Maxwell's equations for a single, 
isolated slit [l^ lT^ll^ . The calculations show a big, up 
to 5 times, resonant transmission enhancement near to 
the Fabry-Perot wavelengths determined by an array pe- 
riod. To clarify the numerical result, we then present an 
intuitively transparent analytical model, which quantita- 
tively explains the resonant enhancement in terms of the 
far-field-zone interference of the evanescent waves pro- 
duced by the independent slits. The model predicts the 
^5-times (resonant) and ~1000-times (nonresonant) en- 
hancements for both the light and matter waves passing 
through a perforated metallic or dielectric screen, inde- 
pendently on the apertures shape. Verification of the 
analytical formulae by comparison with data published 
in the literature supports these predictions. The Wood 
anomalies in transmission spectra of optical gratings, a 
long standing problem in optics follow naturally 

from interference properties of the model. 

Let us first investigate the light transmission by using 
the rigorous model. The model considers an array of M 
independent slits of width 2a and period A in a screen 
of thickness h <^ \. The screen placed in vacuum is illu- 
minated by a normally incident TM-polarized wave with 
wavelength A — 27rc/w = 27r/A:. The magnetic field of the 
wave H{x,y,z,t) — [/(x)exp(— ^(fcz ujt))ey is assumed 
to be time harmonic and constant in the y direction. The 
transmission of the slit array is determined by calculat- 
ing all the light power P(A) radiated into the far-field 
diffraction zone, a;G[— oo, oo] at the distance z 3> A from 
the screen. The total per-slit transmission coefficient, 
which represents the per-slit enhancement in transmis- 
sion achieved by taking a single, isolated slit and placing 
it in an Af-slit array, is then found by using an equation 
TmW = PW/MPi, where Pi is the power radiated by 
a single slit. Figure 1 shows the transmission coefficient 
Tm(A), in the spectral region 500-2000 nm, calculated for 
the array parameters: a = 100 nm, A = 1800 nm, and 
6 = 5 X 10~'^Amaa;- Thc transmitted power was computed 
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FIG. 1: The per-slit transmission Tm{X) of an array of inde- 
pendent slits versus the wavelength for different number M 
of slits. There are three Fabry-Perot resonances at the wave- 
lenghts A„~A/n, n—1, 2 and 3. 



by integrating the total energy flux at the distance z — 
1 mm over the detector region of width Ax = 20 mm. 
The transmission spectra Tm(A) is shown for different 
values of M. We notice that the spectra Tm{^) is peri- 
odically modulated, as a function of wavelength, below 
and above a level defined by the transmission ri(A) = 1 
of one isolated slit. As M is increased from 2 to 10, 
the visibility of the modulation fringes increases approxi- 
mately from 0.2 to 0.7. The transmission Tm exhibits the 
Fabry-Perot like maxima around wavelengths A„ = A/n. 
The spectral peaks increase with increasing the number 
of slits and reach a saturation (TJ^"-^ « 5) in amplitude 
by M = 300, at A w 1800 nm. The peak widths and the 
spectral shifts of the resonances from the Fabry-Perot 
wavelengths decrease with increasing the number M of 
slits. From the data of Fig. 1, one can understand that 
enhancement and suppression in the transmission spectra 
are the natural properties of an ensemble of independent 
subwavelength slits in a thin (6 ^ A) screen. The spectral 
peaks are characterized by asymmetric Fano-like profiles. 
Such modulations in the transmission spectra are known 
as Wood's anomalies. The minima and maxima corre- 
spond to Rayleigh anomalies and Fano resonances, re- 
spectively 'l5]. The weak Wood's anomalies are present 
also in the case of M = 2, a classical Young type two-slit 
system. 

The above-presented data is based on calculation of the 
energy flux by using the electromagnetic field evaluated 
numerically. The transmission enhancement is achieved 
by taking a single, isolated slit and placing it in an ar- 
ray. The interference of the waves diffracted by the inde- 
pendent slits can be considered as a physical mechanism 
responsible for the enhancement. To clarify the numer- 
ical results and gain physical insight into the enhance- 



ment mechanism, we have developed an analytical model, 
which yields simple formulae for the diffracted flelds. For 
the fields diffracted by a narrow (2a -C A, 6 > 0) slit into 
the region \z\ > 2a, it can be shown that the Neerhoff and 
Mur model simplifies to an analytical one. For the mag- 
netic H = (0, iJy,0) and electric E — {Ex,0,Ez) fields 
we found: 

Hy{x, z) = iaDF^{k[x'^ + z^^^), (1) 

E,{x, z) = -az[a;2 + z^]-^/^DFlik[x'' + z^]!/^), (2) 
and 

E,{x, z) = ax[x^ + z^]-'^^DFlik[x^ + z^]^/^), (3) 
where 

D = 4k-^[[exp{ikb){aA - k)]"^ - {aA + kf]'^ (4) 

and 

A = F^(ka) + ^[Fo(ka)Fl{ka) + Fi{ka)F^{ka)]. (5) 

Here, F^, Fq, Fq and Fi are the Hankel and Struve 
functions, respectively. The fields are spatially nonuni- 
form, in contrast to a common opinion that a subwave- 
length aperture diffracts light in all directions uniformly 
16]. The fields produced by an array of M independent 
slits are given by E{x, z) — X]m=i ^mix + mA, z) and 
H{x, z) — X]m=i Hm{x -t- mA, z), where Em and Hm are 
the fields of an m-th beam generated by the respective 
slit. As an example. Fig. 2(a) compares the far-field dis- 
tributions calculated by the analytical formulae (1-5) to 
that obtained by the rigorous model. We notice that the 
distributions are undistinguishable. Thus, the analytical 
model not only supports results of our rigorous model, 
but presents an intuitively transparent explanation of the 
enhancement in terms of the interference of the fields 
produced by the multi-beam source. The array-induced 
decrease of the central beam divergence (Figs. 2(b) and 
2(c)) is relevant to the beaming light 17], and the non- 
diffractive light and matter waves jl8l] . 

The analytical model accurately describes the fields 
E and H , but what is about the transmission coeffi- 
cient? The field power P, which determines the co- 
efficient T^v/, is found by integrating the energy flux 
\S\ ^ \E X H* + E* X H\. Thus, the model accurately 
predicts also the light transmission. We now consider the 
predictions in light of the key observations published in 
the literature for the two fundamental systems of wave 
optics, the single-slit and two-slit systems. The major 
features of a single-slit system are the intraslit resonances 
and the spectral shifts of the resonances from the Fabry- 
Perot wavelengths . In agreement with the predictions 
0, the formula (4) shows that the transmission T = 
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FIG. 2: Electromagnetic fields in the far-field zone, (a) The 
fields Re{E^{x)) [A and D), RB{Hy{x)) {B and E), and 
Re{10Ez{x)) (C and F) calculated for M = 10 and A = 
1600 nm. The curves A, B, and C: rigorous model; curves 
D, E, and F: analytical model, (b) Re{Ex(x)) for M=\: 
analytical model, (c) Re{Ex(x)) for M=5: analytical model. 
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FIG. 3: The per-slit transmission coefficient T (X) versus 
wavelength for the Young type two-slit experiment fl9l| . Solid 
curve: experiment; dashed curve: analytical model. Parame- 
ters: a = 100 nm, A = 4900 nm, and h = 210 nm. 



P/Po = {a/k)\R.e{D)f + \l-m{D)Y exhibits Fabry-Perot 
like maxima around wavelengths A„ = 2h/n, where Pq is 
the power impinging on the slit opening. The enhance- 
ment and spectral shifts are explained by the wavelength 
dependent terms in the denominator of Eq. (4) . The en- 
hancement (T(Ai)«6/™ Hi) is in contrast to the atten- 
uation predicted by the model Q . The Young type two- 
slit configuration is characterized by a sinusoidal mod- 
ulation of the transmission spectra T2(A) 0,1131. The 
modulation period is inversely proportional to the slit 
separation A. The visibility V of the fringes is of order 
0.2, independently of the slit separation. In our model, 
the transmission is given by j[Fl{xi)[iF^{xi)\* + 
Fl{x2)*iFQ{x2)]dx, where xi — x and X2 = x + K. The 
high-frequency modulations with the sideband-frequency 
/s(A) «/i(A) + /2(A,A)~1/A (Figs. 1 and 3) are pro- 
duced as in a classic heterodyne system by mixing two 
waves having different spatial frequencies, /i and /2. 
Although our model ignores the plasmons, its prediction 
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FIG. 4: The per-slit transmission Tm{X) versus wavelength 
for the different values of A and M: (A) A = 100 nm, M = 2; 
(B) A = 500 nm, M = 2; (C) A = 3000 nm, M = 2; (D) A = 
100 nm, M = 5; (E) A = 500 nm, M = 5; (F) A = 3000 nm, 
AI = 5. Parameters: a = 100 nm and b = 10 nm. There are 
two enhancement regimes at A ^ A and A>A. 



for the visibility {V ~ 0.1) of the fringes and the resonant 
wavelengths A„ = A/n compare well with the plasmon- 
assisted Young's type experiment li] (Fig. 3). Some 
difference between the calculated and measured values 
shows that the plasmonless and plasmon-assisted reso- 
nances can compete in certain situations. In the case of 
b > A/2, the resonances at A„ = A/n can be accompanied 
by the intraslit resonances at A„ — 2b /n. 

In order to gain physical insight into the mechanism 
of plasmonless enhancement in a multi-slit (M>2) sys- 
tem, we now consider the dependence of the transmis- 
sion Tm(A) on the slit separation A. We assume that 
the slits are independent also at A ^ 0. According to 
the Van Citter-Zernike coherence theorem, a light source 
(even incoherent) of radius r = M{a + A) produces a 
transversally coherent wave at the distance z<TTRr/X in 
the region of radius R. Thus, in the case of A <C A, 
the collective emission of the ensemble of slits gener- 
ates the coherent electric and magnetic fields, E — 
Sm=i ^■mGxp{i(pm)^MEiexp{iip) and HKMHiexp{iLp). 
Consequently, the maximum power of the emitted light 
scales with the square of the number of slits (beams), 
P ^ AP. Therefore, the transmission {Tm^P/M) grows 
linearly with the number of slits, Tm^M . For a given 
M, the function Tm{}<) monotonically varies with A. 
Such an enhancement regime (first regime) is shown in 
Fig. 4. At appropriate conditions, the transmission can 
reach 1000-times enhancement (M — \z/TiR{a + A)). 
The transmission enhancement in the second regime does 
not require close proximity of the slits. In the case of 
R > Xz/nr (A>A), the beams arrive at the detector with 
different phases ipm- Consequently, the power and trans- 
mission grow slowly with the number of slits (Figs. 1- 
4). According to our model, the transmission Tm ex- 
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hibits the Fabry-Perot like maxima around wavelengths 
An = A/n. The constructive and destructive interfer- 
ence of the beams leads respectively to the enhancement 
and suppression of the transmission amplitudes as in a 
classical heterodyne system. 

The analytical model gives not only intuitively trans- 
parent explanation of the plasmonless transmission en- 
hancement, but in contrast to the previous studies, pre- 
dicts the enhancement for the matter waves. Indeed, in 
our model, the enhancement is based on coherent excita- 
tion of an assemble of slits and constructive interference 
of the diffracted waves. The constructive interference is 
provided not by coupling between the slits, but by a geo- 
metrically well-defined phase relationship between wave 
amplitudes at different lateral locations in the far-field 
zone. Thus, the enhancement mechanism depends nei- 
ther on the nature (light or matter) of the fields -ipm nor 
on material and shape of the apertures. For instance, in 
the first enhancement regime, the fixed phase correlation 
leads to the M-times enhancement of the field ampli- 
tude, V = Z)m=i V"™ exp(i(^m)«A^V'i exp(«(/j), and con- 
sequently to the M-time transmission enhancement. Due 
to Babinet's principle, the model predicts the enhance- 
ment also in the refiection spectra. The destructive in- 
terference of the fields at the detector can lead also to the 
zero transmission. Indeed, the interference of the positive 
(-0) and negative {—ip = ipe'scpiiip + tt)) fields produces 
a field with the zero amplitude and energy. The value 
T=0 is obtained by summing up the respective positive 
and negative energies. Notice that the amplitudes of the 
fields Ipm can rapidly decrease with increasing the dis- 
tances X, y and z. However, due to the enhancement 
and beaming mechanisms (Figs. 1-4), an array produces 
a nonevanescent (propagating) wave ip with low angle 
divergence. Such a behavior is in agreement with the 
Huygens-Fresnel principle, which considers a propagating 
wave as a superposition of secondary spherical waves. 

It is worth noting that the presented model is similar 
in spirit to the dynamical Bloch-waves diffraction model 
|2l| . the Airy-like model based on the Rayleigh field ex- 
pansion [23, and especially to a diffracted evanescent 
wave mode l llll. The difference arises from a fact that 
the models |llLl2llE3 | consider the case of M = 00, while 
predictions of our model strongly depend on the number 
M ^ 00 of slits. In addition, our model deals with in- 
dependent slits, while the models of Refs. [ill EH 
consider the slits electromagnetically coupled via the pe- 
riodic boundary conditions. There is an evident resem- 
blance also between our model and a Dicke superradi- 
ance model of collective emission of an ensemble of 
atoms. A quantum reformulation of our model can help 
us to understand why a quantum entangled state of pho- 
tons is preserved on passage through a hole array [2J|. 
The quantum model will be presented in our next paper. 
Notice, that the surface waves can couple the radiation 
phases of the different slits, so that they get synchro- 



nized, and a collective emission can release the stored 
energy as an enhanced radiation. This kind of enhance- 
ment is of different nature compared to our model. The 
model does not require coupling between the slits and 
does not contain surface waves. 

In conclusion, a rigorous model based on a Green's 
function formalism and an analytical model were pro- 
posed for description of the transmission of light through 
an array of independent slits. Using the analytical model 
we showed a physical mechanism that provides a big 
transmission enhancement of light or matter waves not 
by coupling to the surface waves but by the interference 
of diffracted evanescent waves in the far-field zone. The 
verification of the analytical formulae by comparison with 
data published in the literature supports the predictions. 
The Wood anomalies in transmission spectra of optical 
gratings, a long standing problem in optics, follows nat- 
urally from interference properties of the model. The an- 
alytical formulae can be useful for experimentalists who 
develop nanodevices based on transmission and beaming 
of light or matter by subwavelengths apertures. 

This study was supported by the Hungarian Scientific 
Research Foundation (OTKA, Contract No T046811). 
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